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ABSTRACT 

We report on the detection of an emission line near 1.3 keV, which we as- 
sociate with blue-shifted hydrogen-like sulfur (S XVI), in a 76.8 ksec Chandra 
HETGS spectrum of the afterglow of GRB 020813. The line is detected at 3.3a 
significance. We also find marginal evidence for a line possibly due to hydrogen- 
like silicon (Si XIV) with the same blue-shift. A line from Fe is not detected, 
though a very low significance Ni feature may be present. A thermal model fits 
the data adequately, but a reflection model may provide a better fit. There is 
marginal evidence that the equivalent width of the S XVI line decrease as the 
burst fades. We infer from these results that a supernova likely occurred ^ 2 
months prior to the GRB. We find no discrete or variable spectral features in the 
Chandra HETGS spectrum of the GRB 021004 afterglow. 

Subject headings: gamma rays: bursts — supernovae: general — X-rays: general 
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1. Introduction 

Much evidence has accumulated in recent years connecting long-duration 7-ray bursts 
(GRBs) (Kouveliotou et al. 1993) to supernovae (SNe). The detection of SN 1998bw in asso- 
ciation with GRB 980425 (Galama et al. 1998) and the detection of late-time SN "bumps" in 
several afterglow light curves (Bloom et al. (1999), Reichart (2001) and references therein), 
the observed positional coincidences between star-forming regions and long-duration GRB 
afterglows (e.g. Sahu et al. (1997), Kulkarni et al. (1998), Kulkarni et al. (1999)), and the 
evidence for dust extinction in some GRB afterglows all point toward this association. In 
this context, X-ray spectroscopic observations of early GRB afterglows are tremendously 
interesting, because they potentially allow us to view directly the effects of the GRB on the 
progenitor stellar material and possibly on metals freshly synthesized in the SN. 

Several authors have claimed detection of Fe lines in GRB X-ray afterglows (Piro et al. 
(1999), Piro et al. (2000), Antonelli et al. (2002)). Beyond very weak evidence in Piro et 
al. (2000) for an emission line from H-like S in GRB 991216, Reeves et al. (2002) (hereafter 
R02) are the first to strongly assert evidence for emission from low-Z, multiple-a elements in 
a GRB X-ray afterglow. If valid, this claim has broad implications for GRB emission models 
and it would strongly link GRBs to SNe. The reduction of these data and the statistical 
significance of the results has been called into question by Borozdin & Trudoyubov (2002) and 
Rutledge & Sako (2002), respectively. In a follow-up paper (Reeves et al. 2002b), the XMM- 
Newton ((Jansen et al. 2001)) observers address those concerns. Physical interpretations 
aside, the R02 results are controversial largely due to the small number of observed counts and 
the low spectral resolution of the EPIC-pn instrument, which would blur some fraction of the 
line emission into the continuum. These hinder an accurate measurement of the continuum 
emission and make it difficult to gauge the significance of discrete spectral features. 

We find possible evidence for Si and S emission in the Chandra X-ray Observatory 
(Weisskopf et al. 2002) High Energy Transmission Grating Spectrometer (HETGS) spectrum 
of the GRB 020813 afterglow. We detect a Ni-K line at very low significance, and we find 
a strikingly low upper limit to the Fe-K line equivalent width. The emitting material is 
blue-shifted with velocity ~ 0.1c relative to the GRB host galaxy. These results add weight 
to the detections claimed by R02 and Watson et al. (2002). With the HETGS's high spectral 
resolution, we are able to resolve the line widths, and we are able to accurately separate the 
line and continuum emission. Because the lines we associate with S and Si are detected in 
the full data set, we have ~ 5 times more counts to work with than R02, where the line 
emission is observed during only a short period of time. Moreover, the continuum is an order 
of magnitude stronger than that reported by R02, possibly indicating, in the case of thermal 
emission, that a power-law component is also present. This may alternatively be a sign that 
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the lines are produced under reflection. 

2. Observations 

Using Directors Discretionary Time, two afterglows of GRBs detected and localized by 
the High-Energy Transient Explorer (HETE) satellite (Ricker et al. 2003) were observed with 
the Chandra HETGS, approximately 1 day after each burst, and lasting for approximately 
1 day. Both bursts were long-duration GRBs, with 7— ray durations ^100s. The 7— ray 
peak flux from GRB 020813 (Villasenor et al. 2002) was quite high (~ 10~ 5 erg cm" 2 s _1 ), 
while the peak flux from GRB 021004 (Shirasaki et al. 2002) was approximately two orders 
of magnitude lower. 

Chandra began observing the GRB 020813 afterglow with the HETGS on August 13.990 
(tburst + 21.02 hrs) and continued until August 14.892 (tb urst + 42.67 hrs), with a total 
livetime of 76.8 ksec. The mean count rate from the HEG and MEG gratings and the 
Oth-order image was 0.12 counts/s (0.06 counts/s for the gratings alone). The afterglow 
candidate reported by Fox (2002a) was detected and was observed to fade in brightness with 
time according to a power-law of index a = —1.38 ± 0.06 (xl = 36.4/24), consistent with 
the value reported by Vanderspek et al. (2002) and consistent with measured values in the 
optical (Malesani et al. 2002). Between 0.6 and 6 keV, we measure a time-averaged flux of 
2.2 x 10~ 12 erg cm~ 2 s _1 . 

As first reported by Sako & Harrison (2002), Chandra began observing the afterglow 
of GRB 021004 on October 5.358 {t burst + 20.49 hrs) and continued until October 6.400 
(tburst + 45.50 hrs), with a total livetime of 86.7 ksec. The mean count rate from the HEG 
and MEG gratings and the Oth-order image was 0.04 counts/s (0.02 counts/s in the gratings 
alone). The afterglow candidate reported by Fox (2002b) was detected and was observed to 
fade in brightness with time. The X-ray light curve appears variable. (Similar variability has 
been reported in the optical for this afterglow; see e.g. Halpern et al. (2002)). The temporal 
fade in flux can be described by a power-law with index a = —0.9 ± 0.1 (xl = 83.8/27), 
consistent with the values reported by Sako & Harrison (2002) and Holland et al. (2002). 
The time-averaged flux between 0.6 and 6 keV is measured to be 6.3 x 10~ 13 erg cm~ 2 s -1 . 

The total number of counts detected by Chandra m the case of GRB 020813 (GRB 021004) 
is 6.5 times (2 times) the number of counts detected for GRB 991216 (Piro et al. 2000). 
We perform a parallel spectral/temporal analysis below for the X-ray afterglow spectra of 
GRB 020813 and GRB 021004. 
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3. Spectral Fitting 

We reduce the HETGS spectral data from the standard L2 event lists using IDL and 
custom scripts, described in Marshall et al. (2002). For the Oth-order spectral data, we use 
the CIAO 4 processing exclusively. Each data set is corrected for QE degradation 5 due to 
contamination in the ACIS chips, prior to spectral fitting. We also destreak the S4 chip 6 . 
Spectral fitting and analysis is performed with ISIS 7 . For each observation, we fit the HETGS 
HEG and HETGS MEG lst-order data for the entire observation jointly. The +1 and -1 
orders are combined for each grating. The data are then binned to a S/N > 5 per bin, 
with the bin width restricted to SE/E < 0.1. We define S/N as the background-subtracted 
number of counts divided by the square root of the sum of the signal counts and the variance 
in the background. Bins of the maximal width with S/N < 3 are rejected, due to lack of 
signal. For GRB 020813, this selects an HEG energy range of 0.8 — 7.6 keV and and MEG 
energy range of 0.6 — 5.2 keV. For GRB 021004, this selects an HEG energy range of 1.0 — 4.6 
keV and and MEG energy range of 0.6 — 4.7 keV. We fit each model by minimizing x 2 - Unless 
otherwise noted, all quoted errors are 90% confidence. 



3.1. GRB 020813 

As reported by Vanderspek et al. (2002), we find that an absorbed power-law fit to the 
lst-order spectra requires no absorption column in excess of Galactic value. Therefore, in 
each model that we consider, we include absorption with N H frozen to 7.5 x 10 20 cm" 2 . 

The data are moderately well fit by a power-law (Figure 1, xl = 176.0/154), with 
photon index T = 1.85 ± 0.04 and normalization (5.5 ± 0.2) x 10~ 4 photons/keV/cm 2 /s 
at 1 keV. A thermal bremsstrahlung fit is a much poorer fit to the data (xl — 196.7/154, 
rejectable at 99% confidence), with kT = 4.6 ± 0.5 keV. We find consistent fits using the 
Oth-order data. Of note, a large fraction of the contribution to x 2 i n the case of the HETGS 
power-law fit arises between 1.2 and 1.5 keV. The addition of a Gaussian component in 
this energy range improves the x 2 {xl = 160.5/151). This improvement is at the 99.9% 
confidence level according the the likelihood ratio test (A% 2 = 15.5 for 3 additional degrees 
of freedom). However, this method for determining significance is deprecated (see Protassov 

4 http:/ /cxc. harvard.edu/ciao/ 

5 htt p : / / asc . harvard . edu /cal/ Acis / CaLprods / qeDeg / 

6 ht t p : / / asc. harvard . edu /ciao/ threads /spectraJietgacis / 
7 http://spacc.mit.cdu/CXC/ISIS/ 
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et al. (2002)), and we present Monte Carlo calculations for the significance below. The line 
energy is 1.31 ± 0.01 keV, with a width of 8 ± 4 eV {la). The rest-frame equivalent width is 
67 ± 32 eV, and the flux is (1.6 ± 0.8) x 10~ 14 erg cm -2 s -1 . Given the optically measured 
GRB host galaxy absorption redshift z > 1.254 (Price et al. 2002), a likely association for 
this line is S XVI Ka, blue-shifted by 0.12c from z = 1.254. The x 2 is marginally improved 
(xl = 155.7/150) with the inclusion of an additional Gaussian component of the same width, 
constrained to lie at predicted energy of Si XIV Ka (1.01 keV at z — 0.99 ±0.01). From the 
likelihood ratio test, the confidence for the pair can be estimated as 99.96% (Ax 2 = 20.3 
for 4 additional degrees of freedom). Allowing the widths to vary jointly, the best fit width 
value is 10 ± 9 eV. 

In gauging the significance of these features, we must take into account the number of 
trials performed (i.e. line widths and line energies searched). Figure 2 displays the power-law 
model found above over-plotted on the summed HEG+MEG counts, at three binnings. The 
finest binning is chosen so that each bin between 1.5 and 16.0 A contains one or more counts. 
We rebin twice by a factor of approximately two in order to look for broader features. The 
1.31 keV (9.4 A) feature is seen in both the top and middle plots at 4.0o\ This significance 
describes the Poisson probability that the signal plus background counts could reach a value 
greater than or equal to the observed values, with the Poisson mean fixed to the number of 
counts predicted by the model plus the background. The probability for a chance fluctuation 
of this or greater significance in one or more of the 235 bins shown in Figure 2 is approximately 
1.5%, calculated as P N « 1 — (1 — Pi) N . From Monte Carlo simulations of 10 4 spectra using 
the best fit power-law model above, we find P N = 1.4%. From the simulations, we find that 
the probability for getting two or more 4a fluctuations (as observed) is 0.1%. Thus the line is 
detected at 3.3<r significance. A broad feature corresponding to the 1.01 keV line mentioned 
above is observed in the third panel (and also in the 1st panel) of Figure 2 at ~ 12.2 A, 
with a significance of 2.9<r. We estimate a 58% probability that a 2.9<r or greater fluctuation 
would occur by chance in one or more of the trials. From the Monte Carlo simulations, the 
probability for getting two or more 2.9a fluctuations is 23%. Thus the line is detected at 
very low (1. 2a) significance in a blind search, and it is potentially meaningful only if the line 
location is constrained as in the preceding paragraph. No additional features at or above 
2.9a are present in the HETGS spectrum at the binnings searched. Finally, we derive a 3a 
upper limit of 20 eV for the equivalent widths of any lines unresolved in this search (FWHM 
<C 0.11 A) between 6 and 16 A. This figure grows oc A~ 3 shortward of 6 A. 

We also attempt to fit these data using a collisional ionization equilibrium (CIE) plasma 
model, as has been done previously in the case of several recent GRB afterglows observed 
with XMM-Newton (R02, Watson et al. (2002)). We form the model using the Astrophysical 
Plasma Emission Database (APED) accessible through ISIS, with the redshift fixed at z — 
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0.99. We utilize a turbulence velocity vt, fixed at 2000 km/s, and thermal line profiles to 
produce the best-fit line widths above. We fix the abundance of metals lighter than Mg 
to solar. We tie the Si abundance to the S abundance, and allow Ni to vary freely. All 
other metals are frozen to zero abundance. The CIE plasma fit (Figure 3, xl = 177.8/152) is 
slightly poorer than the pure power-law fit. Allowing the APED metal abundances previously 
set to zero to vary does not improve the fit. The best fit plasma temperature is 9.2 ±1.2 
keV. The best fit S,Si abundance is 2.9 times solar, and the linked abundance is greater than 
solar at 98% confidence (Ax 2 = 5.5 for 1 additional parameter). The best fit Ni abundance 
is 1.8 times solar. There is a low (1.6a) significance detection of a Ni-K line in the MEG 
spectrum. This features appears as two bins modestly over the power-law fit in Figures 1, 
3 near 4 keV. The APED model fits these bins with a blend of K lines from H and He-like 
Ni. The APED model with an abundance of Fe comparable to the Si and S abundances 
would produce a very strong and easily detectable Fe-K feature (Figure 3 Inset). In Table 1, 
we report 1-parameter 90% confidence upper limits to the abundances of several elements, 
with upper limits to the line fluxes and rest-frame equivalent widths of H-like K lines in each 
species. In Figures 4, 5, and 6 we compare the abundances, line luminosities, and rest-frame 
equivalent widths, respectively, to those reported for other GRB X-ray afterglows. 

Since, the abundances of a-particle nuclei produced by core-collapse SNe are expected 
to be approximately the same (see e.g. Rauscher et al. (2002)), relative to solar, the low 
relative upper limits we find for the abundances of Ar and Ca possibly indicate that the 
Si and S lines are not the result of thermal emission from a hot plasma. After radioactive 
decay, we would expect an Fe abundance similar to the low-Z, multiple-a metal abundances. 
This implies that the line fluxes from one or more of the Fe group elements (Fe, Ni, Co) 
should be higher than we observe. We discuss below the possibility that the lines are due 
to reflection. An alternative possibility, suggested in part by the APED model's apparent 
tendency to under- fit the continuum below ~ 1.5 keV (Figure 1), is that we must also include 
a power- law component. It is then possible to use the APED model with a lower temperature 
(~ 2 keV), and this weakens the high energy continuum and the lines from Ar, Ca, and the 
Fe group elements relative to the Si and S lines. We find models with a single abundance 
for the a-particle nuclei and Fe which improve upon the power-law fit but do not appear to 
be as good as the power-law plus two Gaussian fit. Finally, we note that the abundances 
quoted in Table 1 lose meaning if a power-law component is also present. 
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3.2. GRB 021004 

An absorbed power-law fit to the lst-order spectra does not require an absorption col- 
umn in excess of the Galactic value, as reported by Sako & Harrison (2002). Therefore, in 
each model that we consider, we include absorption with Ah frozen to 4.24 x 10 20 cm -2 . 

The data are very well fit by a power-law (Figure 7, xl = 39.5/55), with T = 2.01 ±0.08 
and normalization (1.7±0.1) x 10~ 4 photons/keV/cm 2 /s at 1 keV, consistent with the values 
found by Sako & Harrison (2002) and Holland et al. (2002). A thermal bremsstrahlung fit 
is a somewhat poorer, though still very acceptable, fit (xl — 47.2/55), with kT = 3.2 ± 0.4 
keV. We find consistent fits using the Oth-order data. 

We plot the combined HEG/MEG data in count space with the power-law model over- 
plotted in Figure 8, for the same binnings as in Figure 2. No significant features are observed. 
The largest fluctuation occurs in the middle panel, where one of the ~ 70 bins reaches 2.6a. 
There is approximately an 89% chance that this would happen by chance in one or more of 
the bins in that plot. At the 3a level, this search is sensitive to moderately broad (FWHM ~ 
5X) emission and absorption lines between 6 and 16 A with equivalent widths greater than 50 
eV. Over this same wavelength interval, we derive a 3a upper limit of 35 eV to the equivalent 
widths of any narrow lines (FWHM <C 0.11 A). These figures grow oc A~ 3 shortward of 6 A. 
Considering the GRB redshift of 2.328 determined in the optical by Mirabel et al. (2002), 
several K-lines from S, Ar, Ca, Fe, Co, and Ni would have been detected were their equivalent 
widths larger than these estimated limits. 



4. Temporal Analysis 

As displayed in Figure 9, the afterglow continua, though moderately variable in the 
case of GRB 021004, do not appear to evolve spectrally. To search for time-variable discrete 
features, we divide the summed HEG+MEG data sets into 2, 4, or 8 portions of equal 
numbers of counts. For each time section, we examine the counts data at 3 binnings: 0.11, 
0.21, and 0.42 A ( as in Figures 2, 8). For each GRB afterglow, we fit absorbed power- 
law's to each of these 42 data sets by minimizing the Cash statistic (Cash 1979). We then 
search the bins between 1 and 16 A for deviations from the best fit model, interpreting the 
deviations using Poisson statistics as in Section 3. We set a limiting detection threshold of 
3.7a so that fewer than one chance detection is expected in 15 x 235 bins searched. The 
search yields no emission or absorption line detections in the case of GRB 021004, while for 
GRB 020813 only the Si and S lines (Section 3) are detected. The Si XIV Ka line is detected 
at 3.7a (P/v ~ 53%) in the second quarter (12-28.5 ksec) of data, while the S XVI Ka line 
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(3 detections) peaks at 4.3<r (Pn ~ 6%) in the first quarter (12 ksec) of data. In Figure 10, 
we plot the S XVI Ka line flux (measured with Gaussian fits of fixed width and location) 
and the continuum flux versus time. Using the likelihood ratio test, we find evidence for a 
decrease in equivalent width between the first 28.5 ksec of data and the remainder of the 
observation at 97% confidence (Ax 2 = 4.7 for 1 additional parameter). 

5. Discussion 

The most significant discrete feature we observe in afterglow spectrum of GRB 020813 
is clearly the line we associate with S, which we detect at 3.3<r significance. This is a 
somewhat marginal detection. However, the reality of the line is supported by the presence 
of line counts in various times slices throughout the observation and in each of the four 
independent lst-order spectra (MEG/HEG ±1 orders). For the full observation no detector 
anomalies (e.g. flares, flickering pixels) are observed in the ACIS-S CCD regions yielding the 
HETGS MEG/HEG ±1 order counts from this line. The background is negligible between 
0.6 and 6 keV, contributing only 3% of the total counts. Pileup in the grating spectrum is not 
a concern for this observation due to the very low count rate. Excess counts corresponding 
to the line are detected in each independent spectrum, and the numbers are consistent at 
the 1.5a level. The combined lst-order spectra (MEG/HEG) show a consistent number of 
excess counts at the lcr level. The validity of this line detection is possibly strengthened by 
an association with S XVI Ka and by a marginal detection of a line at the expected energy 
for Si XIV Ka. We have performed a parallel reduction for the GRB 021004 afterglow data, 
and this has yielded no convincing evidence for discrete spectral features. 

Although we find possible evidence for temporal variability in the line emission for 
GRB 020813 (Section 4), the lines are detected in the full 76.8 ksec observation, and the full 
data set is moderately well fit by a plasma in collisional ionization equilibrium (Section 3). 
Figures 4-6 show the line and metallicity parameters we derive in relation to those reported 
for other GRBs in the literature. Our required abundances are only modestly super-solar 
(Section 3), in contrast to those reported for several other GRBs. The line luminosities 
are quite low, and the rest-frame equivalent widths are an order of magnitude smaller than 
those reported for any other GRB. Figure 10 shows the GRB 020813 continuum and S XVI 
line fluxes versus time alongside the data for GRB 011211 from R02. In comparison to 
GRB 011211, it appears that the longevity of the line emission (and also the narrowness 
of the lines) was critical to their detection, whereas the 10 times stronger continuum could 
easily have hidden less persistent emission. 

We can interpret the longer duration and later onset of the lines for GRB 020813 relative 
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to GRB 011211 (Figure 10) as a purely geometrical effect due to observing emitting material 
residing at a larger radius R from the GRB. For radially expanding emitting material heated 
by GRB photons or by the GRB shock, the narrowness of the lines requires that the emission 
comes from only a small solid angle (Q <^ 0.4), as viewed from the GRB, consistent with GRB 
beaming models (Frail et al. 2001). The line duration (tu ne ^ 76.8 ksec) can be attributed 
to time delays across a thin shell, and this leads to the requirement R = ^jff ^ 10 16 cm. 
An estimate R ^ 10 18 cm for the maximum likely R is determined using the smallest solid 
angle found for a GRB jet in Frail et al. (2001). For the reprocessing material to lie this far 
out from the GRB, it must have gotten there prior to the GRB. The line identifications we 
have made (Section 3) require a blue-shift ~ 0.1c, and this high outflow velocity implicates 
expanding SN ejecta rather than a stellar wind. The SN would have gone off 0.1 to 10 years 
prior to the GRB. The need for a precursor SN can be avoided if the GRB host (absorption) 
redshift is taken as a lower limit on the redshift, rather than the actual redshift. In that 
case, an identification of the S XVI and Si XIV lines with Ar XVIII and S XVI [z = 1.53, 
no blue-shift), respectively, could be made. 

In the case of GRB 011211, R02 assert that a shell of thermal plasma at R ~ 10 15 
cm can account for the afterglow spectrum and for the duration of the reported emission 
lines. Taking into account radiation transfer across the shell, Lazzati (2002) argues that 
a larger radius is required (R ^ 10 17 cm), because the large afterglow luminosity implies 
a high density, which implies a cooling time much shorter than the observed line duration. 
Following Lazzati (2002), we find that a thermal model for GRB 020813 also implies R ^ 10 17 
cm, except for the case (discussed below) where the luminosity of the thermal component is 
actually overwhelmed by a (~ 10 times brighter) power-law continuum component. In this 
case, a solution with R ~ 10 16 cm is possible. The density is low (n ~ 10 8 cm -3 ), and this 
implies a shell with thickness dR ~ -R/10 for tt = 0.1. At larger radii, the shell must be 
either very thin or very clumped for the contained mass to be similar to that expected from 
SNe (M cjccta £ 2OM , Woosely & Weaver (1995)). 

We see three reasons, which arise from complexities in our modeling, to be cautious in 
interpreting our results as evidence for a persistent thermal component. First, we observe 
lines with best fit widths an order of magnitude larger than the thermal line widths. We 
account for this in the thermal model by introducing a turbulence velocity of 2000 km/s. It 
is also possible that this broadening is kinematical. Second, the best-fit CIE plasma model 
has marginally dissimilar abundances for the multiple-a elements (Figure 4) and very low 
abundances for the Fe group elements relative to the multiple-a elements (Figure 4, Figure 
3 Inset). This would require an unusual chemical composition for the emitting medium, a 
mechanism beyond the APED model for suppressing emission from metals other than Si and 
S, or the addition of a power-law continuum component to the model as discussed in Section 
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3. Third, the data do not tightly constrain the relative contributions to the continuum 
from the CIE plasma and from this possible power-law component. Indeed, the continuum 
appears non-thermal, as would be expected for the X-ray afterglow of the burst. 

Though these considerations do not rule out a thermal model, reflection models may 
more naturally describe a spectrum with broad lines superposed on a power-law continuum 
(see e.g. Ballantyne et al. (2001), Vietri et al. (2001)). In these models, the line broadening 
is due to Compton scattering. If the electron temperature is similar to the line energy, 
we expect ge/E ~ [kT /m e c 2 ) 1 ^ 2 ~ 0.05 (at 1.3 keV) for a single scattering. Because this 
is already a factor ~ 5 times larger than the observed widths, we require the reflector to 
not be highly ionized. Consulting Figure 2 of Lazzati, Ramirez- Ruiz, & Rees (2002), it 
appears that our observed dissimilarity in line luminosities (Table 1, Figure 5) may further 
reinforce this conclusion. Their reflection models with ionization parameter £ in the range 
10 < £ < 100 can produce Si and S line luminosities which are an order of magnitude greater 
than luminosities of lines from Ar, Ca, and Fe. In the case of GRB 011211, Lazzati et 
al. (2002) argue that the large line equivalent widths rule out nearby reprocessor models. 
This is not the case for GRB 020813, and the line duration may be set by a delayed energy 
injection rather than by the geometry of the reprocessing material (see Rees & Meszaros 
(2000), Meszaros & Rees (2001)). In this case, we would not need a two-step (SN then 
GRB) explosion. 

In conclusion, we detect lines from an over-abundance of light metals characteristically 
produced in massive stars during pre-supernova nucleosynthesis in the X-ray afterglow of 
GRB 020813. The long S XVI line duration and narrow line widths can be explained via the 
geometry of the emitting material, and this would likely require a two-step explosion as in 
the supranova (Vietri & Stella 1998) scenario. The time delay between the SN and GRB is 
^ 2 months in this picture. In contrast, our weak detection of a Ni feature, with no detection 
of Fe (or Co), suggests a very short time delay (<; 1 week) between the SN and the GRB. 
As discussed in Woosely, Zhang & Heger (2002), it is not likely possible to observe a GRB 
through a supernova remnant this young; the Ni feature may indicate that the geometric 
picture is wrong. A thermal model adequately fits the full data set. A bright Fe (or Ni or 
Co) feature is not observed, however, and this is perhaps the strongest reason to favor a 
reflection model, where an Fe line can be quenched by Auger auto-ionization (Lazzati et al. 
2002). 

We thank Harvey Tananbaum for his generous allocation of Director's Discretion Time 
for these observations. This research was supported in part by NASA contract NASW-4690. 
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Fig. 1.— The ±1 order HETGS (HEG/MEG) data for GRB 020813 are moderately well fit 
by an absorbed power-law (Section 3). The Oth-order data (S/N > 5 per bin) are also fit 
well by this model. Here we divide the HEG data by a factor of 5 and the Oth-order data by 
a factor of 20 for ease of viewing. 
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Fig. 2. — The best fit power-law model on top of the combined HETGS HEG+MEG data, 
at 3 binnings, for GRB 020813. The significance (reported in a units) of deviations from the 
model are calculated using Poisson statistics. We associate the residuals near 9.4 A in the 
top and middle plots with H-like S. We associate the high bin near 12.1 A in the top and 
bottom plot with H-like Si. 
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Fig. 3. — The power-law fit to the GRB 020813 ±1 order HETGS data is improved by adding 
two Gaussian's at the engergies of blue-shifted (by 0.12c) S XVI Ka and Si XIV Ka. The 
lines and continuum can be modelled less effectively by a CIE plasma model (bottom plot) 
containing Si, S, and Ni. The inset panel in the bottom plot shows the 2-5 keV portion of 
the CIE plasma model fit with (dotted line) and without (1.8 times solar abundance in) Fe; 
an abundance of Fe similar to the abundances for Si and S is strongly ruled out. The HEG 
data are plotted in grey. The MEG data are plotted in black. 
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Fig. 4. — Metal abundances for the APED plasma model for GRB 020813 (solid circles), 
as in Table 1. Also plotted are the abundances reported in the literature for several other 
GRB X-ray afterglow spectra. The values for GRB 011211 come from R02. The values for 
GRB 010220 and GRB 001025A come from Watson et al. (2002). The value for GRB 991216 
has been taken from Piro et al. (2000). Error bars have been plotted where available. 
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Fig. 5. — Isotropic equivalent line luminosities from Gaussian fits (Table 1) for GRB 020813 
(solid circles, host at z — 1.254) and for several other GRBs in the literature. To determine 
luminosity distances, we use a cosmology with Vl m = 0.3, VL^ = 0.7, and h = 0.65. The 
values for GRB 011211 come from R02, with z = 2.14. The value for GRB 000214 comes 
from Antonelli et al. (2002), with z = 0.47. The value for GRB 991216 has been taken from 
Piro et al. (2000), using z = \. The value for GRB 970828 comes from Yoshida et al. (1999), 
using z = 0.33. The value for GRB 970508 comes from Piro et al. (1999), using z = 0.835. 
Error bars have been plotted where available. 
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Fig. 6. — Rest-frame equivalent widths in eV from Gaussian fits (Table 1) for GRB 020813 
(solid circles, z = 0.99) and for several other GRBs in the literature. The values for 
GRB 011211 come from R02, with z = 1.88. The values for GRB 010220 (z = 1.0) and 
GRB 001025A (z = 0.7) come from Watson et al. (2002). The value for GRB 000214 comes 
from Antonelli et al. (2002), with z = 0.47. The value for GRB 991216 has been taken from 
Piro et al. (2000), using z — 1. Error bars have been plotted where available. 
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Fig. 7.— The ±1 order HETGS (HEG/MEG) data for GRB 021004 are very well fit by an 
absorbed power-law (Section 3). The Oth-order data (S/N > 5 per bin) are also fit well by 
this model. Here we divide the HEG data by a factor of 5 and the Oth-order data by a factor 
of 20 for ease of viewing. 
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Fig. 8. — The best fit power-law model on top of the combined HETGS HEG+MEG data, at 
3 binnings, for GRB 021004. No highly significant deviations from the model are observed. 
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Fig. 9. — To look for spectral variability, the HETGS data for each GRB are divided into 8 
regions of equal counts and fit with absorbed power-laws by minimizing the Cash statistic 
(Cash 1979). The above plots show the power-law indeces T and the (energy) fluxes found 
for each fit, with the best-fit T's derived for the entire data set and the best-fit temporal 
fade found for the count rates overplotted as dashed lines. The fluxes have been divided by 
the mean fluxes for each observation reported in Section 2. Error bars are 90% confidence, 
determined from the 1-parameter confidence intervals. 
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Fig. 10.— The points on the left show the Si XIV Ka line flux for GRB 011211 from R02, 
corrected to the redshift of GRB 020813. The integrated afterglow flux for GRB 012111 
is also plotted and extrapolated. We employ a cosmology with Vl rn = 0.3, Vt^ = 0.7, and 
h = 0.65. At a later time (~ 50 ksec later in the GRB frame), we measure a similar line flux. 
However, the GRB 020813 afterglow flux is approximately an order of magnitude higher that 
that of GRB 011211. On the right, we plot the measured line flux from S XVI Ka, and the 
afterglow flux, versus time since GRB 020813. Error bars and the one upper limit are at la 
confidence. 



-24- 



Element 


Abund. 


Flux 


EW [eV] 


Mg 


4.4 


0.5 


13 


Si 


5.8 


1.9 


53 




(2.9l? : g) 


(1.1 ±0.8) 


(30 ±23) 


S 


4.3 


2.3 


72 




MS) 


(1.6 ±0.7) 


(50 ±22) 


Ar 


1.0 


0.2 


9 


Ca 


1.8 


0.4 


15 


Fe 


0.1 


1.5 


100 


Ni 


4.4 


1.5 


110 




(i.8l?:I) 


1 o +0 - 5 

i - u -1.0 


70+37 
<o_ 73 



Table 1: GRB 020813 line emission upper limits at 90% confidence for APED model elemental 
abundances, H-like Ka line fluxes (xl0~ 14 erg cm" 2 s _1 ), and rest-frame H-like Ka line 
equivalent widths. Values determined from power-law plus Gaussian fits. The observer- 
frame line widths were set via dE/E = v t /c as 0.0067. Line centers are set with z = 0.99. 



